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We have studied cosmological constraints on the number and radii of possible large extra dimensions. 
If such dimensions exist, Kaluza-Klein (KK) modes are copiously produced at high temperatures 
in the early universe, and can potentially lead to unacceptable cosmological effects. We show that 
during reheating, large numbers of KK modes are produced. These modes are not diluted completely 
by the entropy production during reheating because they are produced non-relativistically. This 
means that the modes produced during reheating can easily be the dominant component. For 
instance, for two extra dimensions the bound on their radii from considering only the thermally 
produced KK modes is R < 1.1 x 10~ 4 mm. If the modes produced during reheating are also 
accounted for, the bound is strengthened to R < 2.2 x 10~ 5 mm. This bound is stronger than all 
other known astrophysical or laboratory limits. 

PACS numbers: ll.10.Kk, 98.70.Vc, 12.10.-g 



I. INTRODUCTION 

In the past few years there has been an enormous in- 
terest in the possibility that the presence of large extra 
dimensions can explain the hierarchy problem 0-0], the 
fact that the energy scale for gravitation (the Planck scale 
~ 10 19 GeV) is so much larger than that for the standard 
model (100 GeV). The idea is that the standard model 
fields are located on a 3+1 dimensional brane embed- 
ded in a higher dimensional bulk, where only gravity is 
allowed to propagate. 

This already puts stringent constraints on the size of 
the extra dimensions. Newtons law should definitely hold 
for any scale which has so far been observed. At present 
the best experiments have probed scales down to about 
1mm. Thus, if there are extra dimensions, they can only 
appear at a scale smaller than that. For simplicity we 
make the assumption that the n new dimensions form 
an n-torus of the same radius R n in each direction Q If 
there are such extra dimensions, the Planck scale of the 
full higher dimensional space, Mp jn+ 4, can be related to 
the normal Planck scale, Mp t ±, by use of Gauss' law |Q 

M PA = R n M^ +4 , (1) 

and if R is large then Mp jn+ 4 can be much smaller than 
Mp. If this scenario is to solve the hierarchy prob- 
lem then Mp }Tl+ 4 must be close to the electroweak scale 
(Mp iTl+ 4 < 10 — 100 TeV), otherwise the hierarchy prob- 
lem reappears. This already excludes n = 1, because 
M P . H+A ~ 100 TeV corresponds to R ~ 10 8 cm. How- 
ever, n > 2 is still possible, and particularly for n = 2 



there is the intriguing perspective that the extra dimen- 
sions could be accessible to experiments probing gravity 
at scales smaller than 1 mm. From this point on we use 
M instead of Mp jH+ 4 to simplify notation. 

So far, the strongest constraints come from the obser- 
vation of the neutrino emission of SN1987A @-||. In the 
standard model, a Type II supernova emits energy al- 
most solely in the form of neutrinos. Furthermore, the 
observed neutrino signal fits very well with the theoret- 
ical prediction. If extra dimensions are present, then 
the usual 4D graviton is complemented by a tower of 
Kaluza-Klein states, corresponding to the new available 
phase space in the bulk. Emission of these KK states 
can potentially cool the proto-neutron star too fast to be 
compatible with observations. This has lead to the tight 
bound that R < 0.66/im (M > 31 TeV) for n = 2 and 
R < 0.8 nm (M > 2.75 TeV) for n = 3 §. In fact, an 
even stronger constraint can be obtained from consider- 
ing the contribution to the diffuse gamma background 
from decays of the KK modes produced in cosmological 
supernovae. A conservative estimate yields a bound of 
R < 0.09/xm (M > 84 TeV) for n = 2 and R < 0.19 nm 
(M > 7 TeV) for n = 3 g(J. 

The other obvious place in astrophysics to look for 
these extra dimensions is cosmology p jll]^2]| (see also 
|l3|). In the present paper we go through the possible 
cosmological effects from the presence of large extra di- 
mensions. We solve the Boltzmann equation for the pro- 
duction of KK modes, both during the radiation domi- 
nated epoch and during the reheating phase preceding it. 
We show that unless the maximum temperature reached 
during reheating is very low, the constraints from cos- 
mology are much stronger than the supernova bounds. 



"This assumption has been made in practically all works on 
the subject, however see Ref. M for a different model. 



II. BOLTZMANN EQUATIONS 

The fundamental equation governing the evolution of 
all species in the expanding universe is the Boltzmann 
equation @, L[f] = C[f], where L = df/dt-pHdf/dp 
is the Liouville operator and C is the collision opera- 
tor describing all possible interactions. / is the distri- 
bution function for the given particle species. In the 
present case, there are two terms contributing to the col- 
lision operator: production and decay. There are several 
possible production channels H) : gravi-Compton scatter- 
ing, pair annihilation and bremsstrahlung. In a super- 
nova, nucleon-nucleon bremsstrahlung is by far the dom- 
inant mechanism because of the very high nucleon den- 
sity. However, this is not the case in the early universe, 
the reason being that the early universe is a high entropy 
environment (rj = ns/n^ ~ 10~ 10 ) H]. Therefore NN 
bremsstrahlung is suppressed by a large numerical factor 
~ n 2 N /n 2 . The dominant processes are instead the pair 
annihilation reactions 2j — » KK , vv — > KK, e + e~ — > 
KK Qjll]]. The matrix element for each of these pro- 
cesses is given simply by \M\ 2 — AiS 2 /AMp, where 
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2.4 x 10 18 GeV 



is the reduced Planck mass. 

In order for the standard cosmological equations to ap- 
ply, it is a necessary condition that Pkk *C Pi, where i 
denotes fields on the brane. This means that we can com- 
pletely neglect inverse processes in our treatment. It also 
has the big advantage that we can neglect Pauli blocking 
and stimulated emission factors in the Boltzmann equa- 
tion. In this case we can use the integrated Boltzmann 
equation which is much simpler than the full Boltzmann 
equation JlSj 
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decay, m -> 



(2) 



where m is the mass of the KK state. For the relatively 
low mass modes we look at, the decay lifetime is very 
long PJ. Therefore decays can be completely neglected 
at early times and the production phase can be separated 
from the decay phase. The production equation is then 
given by 

11to 5 T 

h m = -3Hn m + ^pfifm/T), (3) 

where K x (x) is a modified Bessel function of the second 
kind and we have assumed that m e = 0. This assumption 
has very little influence on the results. 



A. Production during the radiation dominated epoch 

The universe enters the radiation dominated epoch at 
some temperature T, which we shall refer to as the re- 
heating temperature, Trh- Production of KK modes dur- 
ing this epoch was studied in detail by Hall and Smith 



pT| , and in this section their results are rederived. The 
present day number density can be found be integrating 
the Boltzmann equation 
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q*K x {q)dq. 
(4) 



This equation applies to the number density for one mode 
with mass to. However, if we are interested in the total 
present day contribution to the mass density from all 
modes, then we need to sum over all modes. This sum 
can be replaced by an integral over dm because the mode 
density is very high M . This integration yields the result 
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where 5 n _i = 2W 2 /T(n/2). In Fig. 1 we show the 
constraints on M from demanding that po < p cr u, for 
the case of n = 2. This result is identical to what was 
found in Ref. @. 
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FIG. 1. The lower bound on M as a function of Trh, when 
only modes produced during the radiation dominated epoch 
are considered. The calculations assume that h = 0.75 and 
n = 2. 

From this it is evident that for Trh ^ 3 MeV the 
bound is tighter than what is found from SN1987A. How- 
ever, it was shown in Refs. [p|JT7| that T RH = 0.7 MeV 
can still be compatible with BBN, so if Trh is sufficiently 
low the overproduction of KK states can be avoided. 
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B. Production during reheating 

In the above treatment it was assumed that the uni- 
verse enters the radiation dominated epoch instanta- 
neously at the reheating temperature. However, this is 
not the case for any physically acceptable scenario. Plau- 
sibly, the universe enters the radiation epoch after some 
reheating by the decay of a massive scalar field (or by 
some other means of entropy production). The only rea- 
sonable alternative is that the radiation dominated epoch 
extended to much higher temperatures (of the order M). 
Here, we look at the "standard" case where reheating oc- 
curs from the decay of the inflaton field (for further dis- 
cussion of inflation in scenarios with large extra dimen- 
sions, see for instance Ref. |l8j and references therein). 

What happens is that the universe starts reheating 
when the inflaton enters the oscillating regime. The im- 
portant parameters are the density, p^.i, of the inflaton 
when reheating begins and the decay rate of the inflaton, 
r^. The Boltzmann equations for this system have been 
solved numerous times (see e.g. Refs. Jl^Jl^]). The result 
is that the temperature of the produced radiation imme- 
diately increases to a maximum value which depends on 
Ptj,^. After this, there is a period of continual entropy 
production during which the universe is matter domi- 
nated by the 4> field and T oc i -1 / 4 (as opposed to the 
case where no entropy is produced, T cx i -1 / 2 ). At the 
time t ~ r^ 1 the inflaton decays rapidly and the uni- 
verse becomes radiation dominated. Using this, it is easy 
to calculate the number of KK modes produced during 
the reheating phase. is directly related to Trh by 
Trh — 0.5y^T\pMp |l7) , but the additional parameter 
p^.i is introduced in the analysis. However, instead of this 
we use the parameter a = Tmax /Trh, where Tmax is 
the maximum temperature reached during reheating (the 
relation between and Tmax is given in Ref. P0f). 
From this, one gets an expression completely equivalent 
to Eq. (|E|). During reheating entropy is produced contin- 
uously. The entropy density is given by s = g*T 3 at all 
times, where g* is the number of relativistic degrees of 
freedom contributing to the entropy jl5| . We assume that 
g» ~ 10.75 at all times. Although this is not the case at 
very high temperatures, the assumption introduces only 
a modest error. Since T oc i^ 1 / 4 and a oc t 2 ! 3 (because 
the universe is matter dominated), the number density 
of a KK-mode (if one ignores production) is n m oc T 8 . 
This yields n m /s oc T 5 during reheating. In order to 
solve the Boltzmann equation during reheating we intro- 
duce the variable X m = n m s _1 (T/T^^) -5 which is con- 
stant during reheating if there is no production. Then 
the Boltzmann equation can be simply written as 
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where we can approximate 

t RH by t RH ~ 1.5g7 1/2 M P T^ [0. The present day 
number density is then given by n mi o = X m ^RHSQ. The 
total density of all KK-modes can be found be integrating 
over dm p"l|], and gives 
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The factor 
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*° comes from the fact that some of 



the produced KK-modes decay before the present. This 
is one of the main new result of the present paper. No- 
tice that Eq. (|§|) is very similar to Eq. (0), except for 
the fact that z n+1 q 3 is changed to z n ~ 6 q 10 . This can 
be understood quite simply. Modes produced during 
reheating are diluted by entropy production relative to 
modes produced during the thermal epoch by the factor 
(T/Trh)~ 5 . Furthermore, the time-interval for produc- 
tion during reheating is different from what it is dur- 
ing the thermal epoch. During reheating dt/dT oc T~ 5 , 
whereas during the thermal epoch dt/dT oc T~ 3 . Alto- 
gether this gives a factor (T/Trh)~ 7 = z~ 7 q 7 difference, 
which is what is found. 

Fig. 2 shows the lower bound on M for n = 2 and 
3, as a function of Trh and a. It is evident from the 
figures that the contours quickly become vertical for in- 
creasing a. This happens because the KK modes pro- 
duced at early times during reheating have been diluted 
away by entropy production so that they do not con- 
tribute significantly at late times. Even so, accounting 
for the modes produced during reheating does strengthen 
the bound on M. For the case of n = 2 the bound is 6.8 
TeV if only the thermal production is included, whereas 
it is 9.7 TeV if reheating is also accounted for (assuming 
Tmax > Trh). For n = 3 these figures are 0.44 TeV 
and 0.66 TeV respectively. 

Notice that this result is different from the result found 
by Giudice et al. |2C|| , that the final abundance depends 
only on Trh , and not on Tmax- The reason is that there 
is a dense spectrum of modes with different masses. The 
abundance of each mode obeys the relations found in 
Ref. PC|| , i.e. that the final abundance for any mode with 
m < few x Tmax is independent of Tmax- However, for 
higher Tmax , many more modes are excited, and the end 
result is that the total production of KK modes increases 
with increasing Tmax ■ 

Nevertheless, if we only consider the bound on M 
from the demanding that po < Pcrit, it is significantly 
weaker than the bound from SN1987a. In the next sec- 
tion we consider the much stronger bound from the dif- 
fuse gamma background. 



3 



8 3 r 




a 3 



\ 

2 3 4 5 

T RH (MeV) 



FIG. 2. The lower bound on M/TeV as a function of T RH 
and a, from demanding that po.thormai + Pq.rh < Pent- The 
upper panel is for n = 2 and the lower for n = 3. The value 
ft = 0.75 has been used. 



III. CONSTRAINTS FROM THE DIFFUSE 
GAMMA BACKGROUND 



Apart from the production mechanisms there is also 
the possibility that the massive KK states decay into par- 
ticles on the brane. The decay rate for different branches 
has been calculated by Han, Lykken and Zhang The 
decay rate into particles on the brane is of the order 
r ~ m 3 /Mp for any kinematically allowed final state. 
With the temperatures discussed here, the decay life- 
time is on most cases of the same order of magnitude 
as the Hubble time. This means that there will be visi- 
ble effects, especially from the decay contribution to the 
diffuse gamma background. It was shown by Hall and 
Smith (see also Ref. that this leads to a very 

stringent constraint on M, even with Trh = 1 MeV. 
The modes produced during reheating have higher mass 
and therefore much higher decay rates. This means that 



even tighter constraints can be put on M if reheating is 
included. 

For n = 2, the contribution to the diffuse gamma back- 
ground from KK decays is given by 
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where (3 = E/T RH and E{z) = exp(-3.3 x 
lQ- 7 z 3 / 2 T RH MeV t (10 10 y)f3 3 ^). As the upper limit 
for the mass integral, 2 max , we take z max = 2.7 x 
W 3 TnH,MeV , because KK modes above this mass de- 
cay before CMBR formation, and therefore do not con- 
tribute to the present diffuse flux. Note that for a = 1 
this equation is equivalent to what is found in Ref. [ pT| . 
Also note that there is the same factor q 7 z~ 7 differ- 
ence between the thermal (the first part in the bracket) 
and the reheating (the second part) contributions as was 
found between Eqs. (|^) and (||). For higher n one ob- 
tains expressions very similar to Eq. (^|). Observation- 
ally, the diffuse gamma background in the MeV range 
has been measured by the EGRET ||Q (30 - 10 4 
MeV) and COMPTEL H,|| (0.8 - 30 MeV) experi- 
ments. The flux measured by EGRET is approximately 
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= 2.3 x 10" 3 (£;/MeV)- 2 07 MeV~ 1 cm^ 2 s~ i ster- 



and the flux measured by COMPTEL is f§ = 6.4 x 
10" 3 (£;/MeV)^ 2 - 3 MeV~ 1 cm- 2 s- 1 ster- 1 . Demanding 
that 4§ „ „ < 4% , translates into the lower bound on 

dti KK — obs 

M shown in Fig. 3 as a function of Trh and a. Note that 
for the relatively low masses we study here, constraints 
from light nuclei abundances J25| are not very important. 

From this figure it is evident that increasing a leads to 
significantly stronger bounds on M. If Tjih — 0.7 MeV, 
then for n = 2 the bound goes from M > 73 TeV at 
a = 1 to M > 167 TeV at a = 1400 (T MAX = 1 GeV). 
For n — 3 the corresponding numbers are M > 3.9 TeV 
at a = 1 to M > 21.7 TeV at a = 1400, a factor of 
5.6 increase. In most reasonable models it is difficult to 
obtain values of Tmax which are smaller than 1 GeV, 
rather Tmax will usually be much higher than 1 GeV. 

Finally, it should be noted that as n is increased, the 
difference between including reheating and only treating 
thermally produced modes increases. In Table I the lower 
bound on M is shown for different values of n and Tmax 
to illustrate this. 

For masses which are low enough that the KK modes 
have not decayed away before the present, the lower 
bound on M as a function of Trh and Tmax ex- 
hibits a quite simple behaviour. For a = 1, M m i n oc 
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whereas for constant Trh and large a, 
This means that for n > 3 the con- 
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FIG. 3. The lower bound on M/TeV as a function of T RH 
and a, from the diffuse gamma background. The upper panel 



is for n = 2 and the lower for n 
has been used. 



3. The value to = 10 'y 



TABLE I. The lower bound on M in TeV for different val- 
ues of n and Tmax- All values are for Trh = 0.7 MeV. 
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Tmax = 0.7 MeV 


50 MeV 


100 MeV 


1 GeV 


2 


73 


161 


165 


167 


3 


3.9 


16.0 


18.6 


21.7 


4 


0.47 


2.96 


3.75 


4.75 


5 


0.10 


0.89 


1.19 


1.55 



radiation dominated epoch was discussed, and the results 
found in Ref. were rederived. 

It was then shown that if reheating before the ther- 
mal epoch is taken into account, the bound on the fun- 
damental Planck scale, M, is strengthened significantly 
and becomes stronger than the supernova bound, albeit 
the theoretical uncertainty is larger. This is the main 
result of the present paper. 

We showed that if Tmax during reheating is 1 GeV, 
then the bound on M is M > 167 TeV for n = 2 and 
M > 21.7 TeV for n = 3. These bound can be translated 
into upper bounds on the radii of the extra dimensions 
using the relation [ pd| 



Rmm — 2 X 10 



31/n-16 



/ !TeV \ 1 



-2/n 



(10) 



From this, we get R < 2.2 x 10 -5 mm for n = 2 and 
R < 2.5 x 10~ 8 mm for n = 3. 

Note that in the present treatment we have assumed 
that the inflaton only decays to matter on the brane. If 
gravitons are also produced at reheating, the bounds are 
tightened. 

We finish by discussing briefly the few possibilities for 
avoiding the very stringent bound obtained above. In 
our analysis we assumed a toroidal geometry for the ex- 
tra dimensions. However, other choices of geometry lead 
to different spectra of KK modes. As shown in Ref. f|, 
compact hyperbolic manifolds can lead to spectra with a 
lightest mode with m > 10 GeV and large energy spac- 
ing. In this case the cosmological bounds disappear com- 
pletely (as does all other astrophysical bounds). The 
second possibility is that there are more branes embed- 
ded in the bulk. In that case, it is possible that the KK 
modes have fast decay channels to light particles on other 
branes. This would mean that the KK modes could have 
decayed into radiation before the present M . 
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tribution from modes produced during reheating keeps 
increasing with increasing Tmax- However, for very high 
masses the modes have decayed away early on so that 
-M„iin reaches a limiting value for high Tmax- 

IV. DISCUSSION 

We have discussed in detail how KK modes are pro- 
duced in the early universe. First, production during the 
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